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Statistics of Wave Functions in Disordered Systems with Applications to 
Coulomb Blockade Peak Spacing 
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Despite considerable work on the energy-level and wavefunction statistics of disordered quantum systems, 
numerical studies of those statistics relevant for electron-electron interactions in mesoscopic systems have been 
lacking. We plug this gap by using a tight-binding model to study a wide variety of statistics for the two- 
dimensional, disordered quantum system in the diffusive regime. Our results are in good agreement with ran- 
dom matrix theory (or its extensions) for simple statistics such as the probability distribution of energy levels 
or spatial correlation of a wavefunction. However, we see substantial disagreement in several statistics which 
involve both integrating over space and different energy levels, indicating that disordered systems are more com- 
plex than previously thought. These are exactly the quantities relevant to electron-electron interaction effects in 
quantum dots; in fact, we apply these results to the Coulomb blockade, where we find altered spacings between 
conductance peaks and wider spin distributions than traditionally expected. 
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I. INTRODUCTION 

The interplay between electron-electron interactions and 
quantum interference has been a theme in condensed mat- 
ter physics for the last two decades. '•^■''■'' The classic con- 
text is interaction effects in disordered systems.'-^ More re- 
cently, mesoscopic systems, such as quantum dots'* or metallic 
nanoparticles,"^ have been intensively investigated; both dif- 
fusive disordered and ballistic chaotic systems have received 
attention. 

A key quantity in studying such effects is the statistics of 
the single-particle wave-functions as one moves from level to 
level, system to system, or in position space. It is well estab- 
lished that weakly disordered quantum systems, as well as bal- 
listic chaotic ones, display universal statistical behavior. ^•^•*2A 
Universal here means that the properties do not depend on the 
microscopic details of the disorder, such as its spatial coiTe- 
lation function. As the behavior is universal, it can be cap- 
tured by relatively simple models. In fact, many properties 
are described well by random matrix theory (RMT); for those 
which involve the spatial behavior of wavefunctions, a sim- 
ple extension of RMT in which eigenstates are described as a 
superposition of random plane waves (RPW) is accurate. ^-^i^ 
In both computational and experimental results, for instance, 
the probability distribution of the spacing between adjacent 
energy levels and the magnitude of the wave function at a 
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?*■ single point closely match RMT predictions 
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Complications arise as disorder increases and the wave func- 
tions become increasingly non-uniform spatially. Such sys- 
tems have been investigated extensively, including localized 
systems'^-"*-''' with wave functions confined to a small area of 
the system. Diffusive systems of intermediate disorder values 
are of the greatest interest to mesoscopic physics; complexi- 
ties due to incipient localization effects have also been studied 
there .^2i^ 

Despite this recent interest, several of the eigenfunction 
statistics most relevant to problems in mesoscopic physics 
have never, as far as we know, been studied. The interac- 
tion contribution to the energy, for instance, involves sums 



over different energy levels of matrix elements of the residual 
(screened) interaction, each of which involves an integration 
over space. Here we particularly study statistics involving 
both different energies and integration over space. We find 
that these statistics deviate strongly from expectations, indi- 
cating that disordered systems are more complex than previ- 
ously thought. We include an analysis of how these results fit 
with previous experimental and theoretical results, including 
calculations from the supersymmetric CT-model.'*~' We then 
apply our statistics to the Coulomb blockade problem. Param- 
eters of the study have been chosen for the greatest relevance 
to the physical conditions of the Coulomb blockade in semi- 
conductor quantum dots. 



We wish to emphasize that the issue here is not the exis- 
tence of a new regime of behavior but rather new character- 
istics in a regime that has been intensively studied for three 
decades. To make this point, it is necessary to demonstrate be- 
yond a shadow of a doubt that our calculations are performed 
in the familiar regime of parameters. Thus, after explaining 
our methodology in Section HI] we establish in Section [111] 
for simple eigenfunction statistics, a general agreement with 
previous results. The comparison with previous work is by 
no means exhaustive but rather serves to confirm the diffusive 
nature of our system and to demonstrate that we get agree- 
ment with RMT/RPW expectations for these well-researched 
statistics. The core of our paper lies in Section II VI where 
we introduce energy correlation statistics that depart markedly 
from cuiTent analytical predictions, thus demonstrating sur- 
prising complexities underlying the disordered model, even 
for parameters that generate agreement with the simple statis- 
tics investigated in Section |lll] The key numerical results 
are in Figs. I12landll3l Section |V] contains a detailed com- 
parison between our numerical results and the existing an- 
alytic results obtained by the supersymmetric cr-model and 
RMT/RPW methods. Finally, we explore the importance of 
these energy correlation statistics in Section IVII where the 
new results are directly applied to the Coulomb blockade. 



II. METHODOLOGY 

All results in this paper were derived within the Anderson 
model, a standard model for describing disordered systemsiii 
The Anderson model employs a discrete lattice geometry, and 
combines a "hopping" Hamiltonian with a set of uncorrelated 
on-site energies: 

H = Y,\^)e,(^\~VY,i\^){^ + l\ + \^-l){^\). (1) 

i i 

In more than one dimension, the hopping component includes 
transitions to all nearest neighbors. For simplicity, we take V, 
the transition amplitude, and a, the lattice spacing, both equal 
to 1. The ei's are uncorrelated from site to site, and ei/V 
are drawn from a uniform random distribution of width W, 
which measures the disorder strength of the system. To treat 
a magnetic field, the appropriate phase factor can be added to 
Vr^^ Our study concentrates on a two-dimensional rectan- 
gle of size 164 x 264 with hard-wall boundaries. The large 
size allows us to choose parameters such that the disorder is 
weak but the mean free path is less than the system size; hard 
wall boundary conditions are more appropriate for realistic 
quantum dots than other simple possibilities; and the asym- 
metry breaks pseudo-degeneracies in the eigenband. Some 
characteristics of other geometries, such as smaller systems 
or periodic boundary conditions (corresponding to a torus), 
are touched on where appropriate. 

Three main variables can be adjusted to control the physical 
regime: (i) the mean energy around which to draw statistics, 
(ii) the disorder strength, and (Hi) the strength of any applied 
magnetic field. 

First, each fully diagonalized matrix in this study would 
produce 164 x 264 eigenfunctions. Because of the large 
size, computing constraints made it reasonable to analyze only 
eigenfunctions within a narrow band of eigenenergies, here 
expressed in terms of a "filling ratio." This ratio indicates 
the position of the chosen band's central eigenfunction in the 
full eigenenergy band. The presented statistics are averaged 
over both the band's many eigenfunctions and separate dis- 
order realizations. We analyze energy bands centered at 1/25 
and 1/100 filling, with band width about 1/100. The very low 
energy in the second case was chosen to match physical quan- 
tum dots, which typically contain only hundreds of electrons. 
As an added benefit, low energies result in a larger wavelength 
{kpa — 0.72 and 0.35, so that Xp ~ 9a and 18a, respectively), 
compensating for the somewhat arbitrary nature of a discrete 
geometry. 

Second, the disorder strength of the system, measured as 
W , has a profound effect on the system's behavior Our results 
will first demonstrate trends as W is varied, and then focus on 
specific W values. Very small W produces a semi-ballistic 
system, whereas very large W produces localized eigenfunc- 
tions. We are most interested in intermediate disorder, which 
yields a diffusive system. 

Finally, we are also interested in the effect of a small mag- 
netic field, enough to break time reversal invariance but not 
enough to cause well-defined Landau levels, for instance. To 
add a magnetic field B' perpendicular to the system plane, we 



adjust the hopping amplitude in the Hamiltonian by using the 
Peirels substitution:^?^ 



%j}{r) 



^z2^/A'd//v.o^(^)^ A' = B'yx. 



(2) 



This implies that one changes the hopping terms in the x di- 
rection according to 



- 1 



_^±i2TTB'ya/ipo 



^±iyaB/A 



(3) 



where a is the lattice constant, A is the area of the system 
(A/a^ is the number of lattice sites), and B is the magnetic 
field in units of magnetic flux quanta through 27r times the 
area of the system. 

We analyze below both the B — case and a system with 
sufficient magnetic field to break time-reversal invariance. By 
analyzing trends as a function of B, we determined that B — 6 
is sufficient. The results in the presence of a magnetic field 
are applied to the Coulomb blockade problem. 



III. SIMPLE STATISTICS : AGREEMENT WITH 
RMT/RPW 

Before delving into relatively unexplored statistics, we 
want to first determine the parameters corresponding to the 
diffusive regime and confirm that, for simple eigenfunction 
statistics, our computations match previous analytical and nu- 
merical work. We do not carry out a comprehensive com- 
parison with past results but rather present enough to con- 
vincingly show that RMT augmented by RPW or perturbative 
techniques accounts for these properties, all as a prelude to 
the striking disagreement presented in the next Section. For 
simplicity, we start with B = Q and study only eigenfunctions 
in the low-energy range most relevant to mesoscopic physics. 

We begin by establishing the disorder strength W corre- 
sponding to diffusivity. To do this, we consider the inverse 
participation ratio (IPR): 



IPR = A I iVWI^d^r. 



(4) 



A is the area of the system, and we understand the inte- 
gral to be taken as a discrete sum over lattice sites. The 
IPR is closely linked to the degree of localization within the 
system,^"* and is inversely proportional to the volume in which 
the wavefunction is confined. Large IPR values thus corre- 
spond to strongly localized states. It can be derived from 
basic RMT that the IPR for a time-reversal invariant system 
should have a universal mean value of 3.0, with small varia- 
tions from this mean.^^-^^ Chaotic systems, and weakly disor- 
dered experimental systems, do show a mean of 3.0, with a 
nearly symmetrical distribution around this value. Iii2i2£ The 
mean IPR obtained by averaging over disorder realizations has 
been extensively studied numerically for both two- and three- 
dimensional systems. ^^i^"'-^''-^^ The expectation is that the IPR 
should remain roughly constant for disorder values in the dif- 
fusive regime, and then should rise sharply for greater disorder 
values in the localized regime. 
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FIG. 1: The mean inverse participation ratio (IPR) as a function of 
disorder for two energies. The mean IPR increases with disorder; the 
plateau in the middle of each curve, most noticeable for 1/25 filling, 
corresponds to the diffusive regime. RMT predicts a universal value 
of 3.0 (dotted line). The system is a 164 x 264 rectangle, at _B = 0, 
with 5 disorder realizations and about 400 different states used in 
each case. 



Figure[2shows our results for the variation of (IPR) against 
disorder strength for both 1/25 and 1/100 filling. Note three 
important features: IPR is consistently higher for the lower 
energy band, IPR rises with disorder strength, and the pat- 
tern of this trend suggests three separate regimes of behav- 
ior. Between W = 0.1 and 0.3, we see a saturating increase 
in both graphs; between W = 0.3 and 0.5, we see a linear 
trend, especially in the 1/25 case; and for higher W values, 
we see rapidly increasing behavior. It is reasonable to ex- 
pect these three statistically distinct regimes to correspond, at 
least roughly, to the three physically relevant regimes: semi- 
ballistic, diffusive, and localized. Figure [^thus suggests we 
can choose W = 0.2, 0.5 and 0.8 or higher as representative 
values of these three regimes in the 1/25 filling case. 

For further verification that these disorder values corre- 
spond to the diffusive regime, we find the three characteris- 
tic length scales for our system: the electron wavelength, the 
mean free path I, and the localization length, iioc- By com- 
paring the size of the system, L, to these three lengths, the 
various regimes are defined. In the diffusive regime one ex- 
pects 



\f ^l^L<^L 



loc- 



(5) 



In the semi-ballistic regime, the system size is smaller than the 
mean free path, while in the localized regime, the localization 
length is smaller than the system size. 

According to calculations using the Born approximationi^ 
at 1/25 filling, l/a = 2A/W'^, and at 1/100 filling, (./a = 
11.5 /W^. For L, we take the small side of the rectangle, L = 
164a. On the other hand, it is known that Lioc = cNchi, 
where c is a constant found to be about 2. A^'ch, the number 
of transverse channels in the system, is based solely on the 
energy of the system, and is given by iVch = kpL/w. At 1/25 
filling, TVch = 37; at 1/100 filling, TVch = 18. We can thus 
calculate, for the two energy levels, the required W values for 
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FIG. 2: (Color online) Probability density of the spacing between 
neighboring energy levels. The calculated distribution (dashed) 
matches the Wigner surmise (solid) obtained from the orthogonal 
ensemble of RMT. Both the mean and integral are normalized to 1. 
(Filling is 1/25 and B = 0.) 



a diffusive regime: 

1/25 Filling : 

1/100 Fining 



24/1^2 ^ iQ^ ^ 1800/1^2 

=^ 0.38 < P^ < 3.3 (6) 

11.5/W^ < 164 < 420/VF^ 

=^ 0.26 < VT <C1.6 (7) 



This confirms the choice of W = 0.5 for 1/25 filling as be- 
longing in the diffusive regime (Xp = 9a, £ = 96a, Ljoc = 
7200a), and prompts the choice ofW^ 0.35 for 1/100 filHng 
( Af = 18a, ^ = 94a, Li^c = 3400a). 

We can now begin to look at some of the basic statistics 
of our diffusive system and compare them to RMT prediction 
and experiment. For clarity, we will present only the 1/25 fill- 
ing results in the zero magnetic field case to establish general 
agreement, and then proceed to the lower energy when we add 
a magnetic field and begin to concentrate on the energy corre- 
lation statistics relevant to the Coulomb blockade problem. 

One of the simplest and most comprehensively studied 
statistics is the distribution of the spacing between adjacent 
energy levels, P{s)J A basic result of RMT is that P{s) is 
very well approximated by the classic Wigner surmise, given 
in the absence of a magnetic field by 



P{s) 



-TiS^/4 



(8) 



Figure 12] compares this prediction with our calculated energy 
spacing distribution. We see excellent agreement on all parts 
of the graph; the match is equally good for 1/100 filling, not 
shown here. 

Another well-studied statistic is the magnitude of the wave 
function at a single point. In the absence of a magnetic field, it 
is expected to follow the classic Porter- Thomas distribution^— 



Pit) 



1 



= -t/2 + _ 



t^mrrA. 



(9) 



Both ballistic chaotic and weakly disordered systems show 
this behavioriMMmiiill^ The predicted IPR value of 3.0 
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FIG. 3: (Color online) Probability distribution of A\ip{r)\^ in three 
different regimes: semi-ballistic at W — 0.2 (squares), diffusive at 
W = 0.5 (circles), and localized at V[^= 1.2 (pluses). In the diffusive 
case, we see excellent agreement with the Porter-Thomas distribu- 
tion predicted by RMT. The integral of each curve is normalized to 
1. Filling is 1/25 and B = 0. (Because of the hard- wall boundary 
condition, the probabilities were sampled on the inner 3/4 of the 
rectangle.) 



is derived by taking the appropriate moment of this Porter- 
Thomas resuh. Note the general RMT prediction of approxi- 
mately uniformly extended wave functions (in position space) 
such that large wave function amplitudes are exponentially 
rare. 

Figure |5] displays our results for the probability distribu- 
tion of |'(/'('')P in the three different regimes. We see excel- 
lent agreement across the graph for the diffusive case, another 
indication that this system is truly diffusive. As we can see 
from the figure, the higher the disorder value W, the greater 
the prevalence of both very large and very small |V'(^)P val- 
ues, indicating the system is becoming more localized and less 
uniform. 

To transition to system-wide eigenfunction characteristics, 
we first look at the spatial correlation statistic, the correlation 
between sites separated by a certain distance: 

V; Correlation = A{ip{r)ip{r + r')) (10) 

V-^ Correlation = ^^ (V/(r)V'^(r + r')) (H) 

Notice that the IPR is the value of the second correlation func- 
tion at r' — 0. RMT gives an overly simple prediction for 
these correlations: All off-diagonal terms are equivalent and 
have a value consistent with normalization of the wave func- 
tion. A much more useful prediction can be obtained from 
Berry's idea that a wavefunction of a chaotic system can be de- 
scribed as a random superposition of plane waves of the same 
wavenumber but different propagating direction. ''^•^^ Random 
plane wave (RPW) modeling gives predictions for both cor- 
relation functions abova^^i^^ which agree with the perturba- 
tion theory results for the diffusive regime: ^^— 
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FIG. 4: Spatial correlation of the wave function in the diffusive 
regime. The data for {tp{r)tp{r + r')) is close to the random plane 
wave result, especially for small r' . Filling is 1/25, W = 0.5, and 
5 = 0. To avoid boundary interference, we used a 30 x 30 section 
of the rectangle centered 1/4 of the side length from each boundary. 
Correlations were measured in the x-direction from points within this 
section (y-direction correlations are identical). 
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FIG. 5: Spatial correlation of the square of the wave function in 
the diffusive regime. Calculated results for {■ip'^{r)ip'^{r + r'j) are 
in reasonable agreement with the random plane wave result, espe- 
cially for r' < 20. RMT predicts that the correlation function should 
rapidly approach 1. See Figure|4|caption for parameters. 



The t/) correlation should thus approach 0, and the i/i^ corre- 
lation should approach 1 for large r' . Figures|4]and|5ldisplay 
our results for the diffusive regime, and show good agreement 
with the RPW predictions, especially in the tp correlation. At 
large distances, the correlations do appear to converge to the 
expected values, although we note visible discrepancies in the 
■0^ case. 

The final characteristic we wish to consider before moving 
to statistics involving different eigenfunctions is the IPR dis- 
tribution. Previously, we touched on the mean IPR values for 
different parameters, but Figure|6]displays a histogram of IPR 
values for specific system realizations. As disorder increases 
across the three regimes, we see four important effects: the 
distribution gets wider and increasingly asymmetrical, and 
both the median and mode IPR values increase. In the dif- 
fusive regime, the distribution is Gaussian-like near its max- 
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FIG. 6: Distribution of the inverse participation ratio (IPR) for the 
three different regimes (semi-ballistic alW = 0.2, diffusive atW = 
0.5, and localized at W = 1.2). At low disorder, the RMT/RPW 
prediction of a thin, Gaussian distribution is met. With increasing 
disorder, and even in the diffusive regime, the distributions become 
clearly asymmetric and increasingly wide. The integral under each 
curve is normalized to 1. (1/25 filling with B = 0.) 
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FIG. 7: The mean of the Mij as a function of disorder strength. 
When averaged over all wave functions in our energy window (cir- 
cles), the mean rises roughly quadratically with increasing disorder 
strength from 1, the RMT value. In contrast, for consecutive energy 
levels (i — j — 1, squares), the mean is suppressed below 1 at weak 
disorder and shows evidence of the three regimes. (Five disorder 
realizations with B = 0.) 



imum (as it should be given the smallness of the cumulants 
beyond the second one^Ml ), but the tails are clearly asym- 
metrical. In the large IPR tail, we find that the data follow an 
exponential distribution with a decay rate of ~ 8.0 (obtained 
from fitting values larger than 3.4). Such an exponential de- 
cay is expected from calculations using the supersymmetric 
CT-model.^'-^^"*^ For our parameters, the predicted decay rate 
is ~ 9.4 [see e.g. Eq. (3.92) in Ref . 12111 . Considering that in 
our case £ < i while the cr-model is valid for £ <^ L, we 
find the good agreement between the cr-model result and our 
observation above to be another demonstration of the univer- 
sality present in these systems. We will return in Section fVl 
below to make a detailed comparison with the analytic results 
for the magnitude and variance of the IPR in the context of the 
correlation between different wave functions. 

We have now considered five separate statistics of vary- 
ing complexity in the diffusive regime. Energy level spac- 
ing, probability distribution at a point, and the two spatial 
correlation statistics all are in excellent agreement with RMT 
and RPW predictions. We note, however, that the latter three 
statistics depend on individual sites alone and not system- 
wide, or global, characteristics. IPR, a global statistic, seems 
at this level to be within the general framework of expecta- 
tions. 

The establishment of this agreement serves two ends: justi- 
fication that we are dealing with an authentic diffusive sys- 
tem of most relevance to mesoscopic physics, and confir- 
mation that we are getting results for simple statistics that 
match previous work as well as analytical predictions. The 
latter demonstrates explicitly that the unexpected and com- 
plex coiTelations to be encountered in the next Section are 
perfectly consistent with the simple and well-analyzed single- 
eigenfunction statistics familiar in disordered models, and so 
are likely to have wide application in disordered quantum 
physics. 



IV. CORRELATION OF DIFFERENT WAVE FUNCTIONS: 
DISAGREEMENT 

The spatial similarity of different eigenfunctions is cru- 
cially linked to electron-electron interactions and so to the 
statistics of quantum dots in the Coulomb blockade regime. 
AH data presented in this section are based on the quantity 
Mij , defined as 



M,,^A \Mrr Hjirr d'r. 



(14) 



where i and j label different eigenfunctions. Afy thus mea- 
sures the system-wide, spatial correspondence between two 
eigenfunctions in a specific disorder realization. Note that Ma 
is the IPR, discussed in the previous section. 

RMT predicts no coiTelation between different eigenfunc- 
tions in the N ^ oo limit. In this case of uncorrected unifor- 
mity, and within the Weyl approximation 
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IM, 



(15) 



mean Mij should equal 1.0 for all i ^ j, the statistics of the 
Alij should be independent of i or j, and there should be no 
correlation between different Mij . This section demonstrates 
that these basic predictions of RMT and its simple extensions 
(such as RPW) are not met in our diffusive system. 



A. B = 0: Time-reversal invariant 

We begin by studying Mij in the simpler zero-field case. 
Figure0shows the trend in Alij against disorder strength W, 
and also compares the average (on both i and j, in addition 
to disorder) for all i j^ j in the band with the average (on j 
and disorder) for consecutive eigenfunctions. In all cases, the 



mean Afy appears to converge to 1 for very low W values, but 
are consistently higher elsewhere. As in Figure^which shows 
IPR as a function of W, Mij rises with increasing disorder, 
and does so rapidly for consecutive eigenfunctions at 1/100 
filling. Again mirroring the IPR dependence, Mij is higher 
at the lower energy and rises sluggishly at 1/25 filling. Note 
that for both energy levels the points at which the plots for 
consecutive and total Mij cross are near the chosen values for 
diffusivity: W = 0.35 for 1/100 filling and W = 0.5 for 1/25 
filling. Finally, we note that calculations for a torus (periodic 
boundaries in each direction) yielded a mean Mij significantly 
closer to 1.0 for most parameters. 

One remarkable feature of Fig.0is that the mean (over dis- 
order and j) of Mij for consecutive eigenvalues i = j — I is 
substantially different from the result obtained by further av- 
eraging this quantity over the index i, implying some degree 
of disorder and energy dependent correlations between the 
eigenfunction. To investigate this further, Fig.|8]plots mean 
Mij against j — i, indicating the correlation's dependence on 
the nearness in energy of the two eigenfunctions. We guess 
from Figure |8] that the points of intersection in Figure rep- 
resent the disorder levels at which Mij is approximately con- 
stant against the difference in i and j, which we thus presume 
to be a characteristic of diffusive systems. In comparison, the 
semi-ballistic system displays a slightly negative correlation 
for close energies and positive correlation for for distant ener- 
gies. The localized system has a stronger positive correlation 
for close energies — states which are nearby in energy tend to 
overlap in space — which decreases as the energy difference 
increases. Note that this is exactly the opposite of the well- 
known trend in the strongly localized case: when L\oc ^ L, 
states that are close in energy tend to occupy different parts of 
the sample. '^-'^ It is quite surprising that modest changes in 
W could so drastically change how energy correlation oper- 
ates in a system, going from increasing with energy difference 
to staying constant to decreasing and presumably back to in- 
creasing in the strong localization limit. 

We cannot, at this time, claim that we understand the origin 
of the energy correlation between the wave functions leading 
to the behavior observed for the A/^ 's in Figure|S] In analytic 
treatments, the energy scale known as the Thouless energy 
plays an important role: for diffusive systems, ii^Th = 'hD/A 
is the energy scale related to the time tn ~ A/D needed to 
diffuse across the whole system. In both supersymmetric a- 
model and RMT/RPW approaches, the expectation is that the 
mean Mij would be independent of i~ j for energy differ- 
ences less than _ETh and then approach 1 rapidly as a power 
law, oc l/{i — iY [see e.g. Eq. (3.84) in Ref. 12111 . This is 
clearly not the case in our data! We would like however to 
mention that, from a qualitative point of view, the observed 
energy correlation would be compatible with the existence of 
a relatively small number of localized resonances. 

What we mean by a localized resonance is, in a very 
schematic way, what would result from the following picture. 
Assume one can define an approximation H^ of the Anderson 
Hamiltonian Eq. Q such that the eigenstates of iJp can be 
divided in two classes: a vast majority of delocalized states 
■0", and a smaller number of very localized states Lpi such that 
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FIG. 8: (Color online) The dependence of mean Mij on the spacing 
between the two states for the three different regimes (semi-ballistic 
at W = 0.2, diffusive at W^ = 0.5, and localized at W = 0.8). The 
point at which j — i = 0, corresponding to the IPR, is omitted for 
clarity. The mean increases with W , and the effect is magnified for 
close eigenfunctions. Note that the curve for the diffusive system is 
nearly flat. (1/25 filling with B = 0.) 



IPR<p, ^ 1. In fact, states which are in some sense "anoma- 
lously localized" are known to exist from supersymmetric cr- 
model investigations.^"-^' Let us furthermore assume that the 
perturbation V = H — Hq couples the ipi to the ipf with ma- 
trix elements whose typical magnitude v is large compared to 
the mean level spacing A, but small enough in terms of the 
energy spacing of the ifi that these latter remain essentially 
decoupled. 

In such a circumstance, we can model the eigenstates ipi 
of the full Hamiltonian H near the energy e; of (pi using a 
resonant level model, implying that 



ipi = a^ipi " "^ + fr{ei - ei)rittpi 



(16) 



Here, ^f°^°^ is a delocalized wavefunction orthogonal to tpi 
(and not necessarily close to tpf), rji is a fluctuating quantity of 
rms one, and the smooth positive function /r(e — e;) of width 
r ^ t; describes the envelope of the resonance. The resonance 
contains approximatively A^ios ~ F/A levels, and normaliza- 
tion imposes that within the resonance (i.e. for ei — ei < T), 
/p ~ 1/Ni-cs- If A^ros is large enough, the normalization factor 
fli is not very far from one, and we shall drop it from now on. 
If Eq. ( I16t is a good model for the eigenstates of H near the 
resonance, we see that, in this region of energy, the envelop of 
the Alij should be given by (neglecting terms of order 1/Aios) 



(A%) ~ 1 + /2(e, - e,)f^{e, - eOlPR. 



Vl 



(17) 



If IPRip, is not negligible compared to N^^^, such an expres- 
sion provides a mechanism for increasing the mean value of 
Mij, and furthermore would explain that this enhancement is 
larger if i and j are close in energy since this increases the 
probability that they belong to the same resonance. 

In the case of a torus (periodic boundary conditions), we 
found that the enhancement of the Afy's is smaller. We in- 
terpret this as implying that the localized resonances are pref- 
erentially created near the hard wall boundaries of the sys- 
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FIG. 9: (Color online) Energy level spacing distribution in a mag- 
netic field large enough to break time-reversal symmetry (6 </3o 
through 2ti:A). The curves at three different energies, all in the 
diffusive regime, match the corresponding Wigner surmise of RMT 
(solid). Both the mean and integral are normalized to 1. 



tem; this is actually the region where we see the first localized 
states appear as disorder is further increased. 

To finish this subsection, we comment briefly on two points. 
First, we mention that one statistic that may help explain some 
features of A/^ is the correlation of the wave-function with the 
disorder configuration of the system, defined as 



Q -2.(^.(0-^)^0, 



(18) 



where e(r) is the on-site disorder at r and A is the area of the 
system. Our computational results demonstrate that this dis- 
order correlation, negative (for filling smaller than one half) 
for all W values, is proportional to the square of the disor- 
der strength and, in fact, matches the value derived from per- 
turbation theory. Since individual wave-functions in a spe- 
cific system are correlated to the same disorder configuration, 
they will be correlated to one another. It should be stressed, 
however, that the amount of correlation thus induced is much 
smaller than that seen in the A/^ data. 

We finally note in passing that several statistics show peri- 
odic structure as a function of energy. The mean IPR and dis- 
order correlation, for example, show small oscillatory effects 
when plotted against eigenfunction number. As the trends are 
more prominent in the torus or smaller systems, they are most 
likely due to periodic orbit effects. However, persistence of 
this odd behavior when the smallest period is about four times 
the mean free path, suggests some extra relationship may in- 
fluence how wave-functions at specific energies interact. 



B. B 7^ 0: Broken time-reversal symmetry 

Because we want to eventually study interaction effects in 
the Coulomb blockade regime, and for these latter the discus- 
sion of the zero magnetic field case is made more complicated 
by the partial screening of the Cooper channel,-'*^ we shall 



now consider a situation where a magnetic field is applied, 
effectively suppressing the time reversal invariance of the sys- 
tem. For the Coulomb blockade, we should like to choose 
a system in the diffusive regime containing a few hundred 
electron. These considerations lead us to pick the following 
parameters: 1/100 filling corresponding to about level 400; 
W — 0.35; and i? = 6, the field strength at which the statis- 
tics indicating broken time-invariance appear to level off as a 
function of B. We also wish to supplement this system with 
data for an intermediate energy range. Thus, three such sys- 
tems will be investigated: Our favored system at 1/100 filling 
and 14^ = 0.35; a system at 1/25 filling and 1^ = 0.5, for com- 
parison with the 5 = case; and an intermediate system at 
1/75 filling and 1^ = 0.5. 

Our first concern is to establish that the three systems with 
non-zero magnetic field are diffusive. Figure|9]plots each sys- 
tem's energy level spacing distribution along with the relevant 
Wigner surmise (GUE for broken time-reversal invariance). 
We see excellent agreement in all three cases, an indication 
of diffusivity. In addition, note that the discussion of length- 
scales - mean free path, system size, and localization length - 
in Sectionirniholds unchanged for the weak field considered, 
and so suggests diffusivity in each of our three cases. Further- 
more, simple statistics like IPR and probability distribution 
(not shown) match RMT-based predictions reasonably well. 
A final confirmation is assurance that the cyclotron radius, 



Re 



mvp 
eB' 



kpA 



(19) 



is much larger than the system. Indeed, finding kp from the 
dispersion relation yields at S = 6 



Re (1/25 Fining) = 830 

Re (1/75 Filling) = 490 

i?c (1/100 Filling) = 400 



(20) 



all considerably larger than the system size (264) or the mean 
free path. We can thus be reasonably assured of diffusivity in 
these three systems. 

Figure^Jdisplays Mij data against closeness of energy for 
our favored 1/100 filling system: the mean and median are 
compared as well as the standard deviation. All three plots 
have a characteristic shape: a rapid rise as i—j decreases from 
about 50 to 3 followed by a tiny dip from about 3 to 1 . As 
the pattern is closest in form to the slightly localized B — 
case, it is possible that the application of a magnetic field may 
be strengthening prelocalization effects. However, this basic 
form holds at non-zero magnetic field even for weak disor- 
der (not shown), discrediting the notion that we are seeing 
genuine localization at these parameters. Note that the mean 
Mij is consistently higher than the median, and significantly 
larger for close energies. Indeed, the distribution of A/^ is 
strongly asymmetrical for close i, j, and increasingly Gaus- 
sian for more distant eigenfunctions. As for the B = case, 
the observed energy dependence in the mean and distribution 
of Mij are neither seen in nor expected from analytic results 
to date. 
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FIG. 10: (Color online) Mean, median, and standard deviation of 

Mij as a function of the spacing between the two states. The scale for FIG. 1 1 : Covariance of Mhj and Mij (distinct h, i, j) in three dif- 



the standard deviation is on the right, and the point j — iis omitted for 
clarity. The similarity in form of the mean and the standard deviation, 
as well as the increasing disparity between mean and median with 
closeness in energy, are striking. (1/100 filling, VF = 0.35, B = 6.) 



fusive cases as a function of the spacing between i and h, averaged 
over all j in the energy window. For clarity, the first ten points in 
each plot are shown intact, whereas the remaining points are aver- 
aged in groups of ten. Note the qualitative difference at 1/25 filling 
between the zero-field and B = 6 cases. 



C. Interaction Terms 



The Mij enter directly into interaction terms for the 
Coulomb blockade peak spacing through Fj , defined as 



FAn] = 



E 



M. 



Mij 



(21) 



where the Mij average (over disorder) subtracted from each 
Mij is calculated from the same system parameters and eigen- 
functions. Note that Mij depends on the eigenstate indices i 
and j. The disorder average of Fj is thus 0, by definition. 
The specific term included in the peak spacing calculation is 
{Fj+i[n] — Fj[n]), whereby a larger magnitude will lead to 
wider peak spacing and spin distributions. We are thus inter- 

Fj[n]), given by 



ested in the root-mean-square of (-Fj+i [n] 

F'^ln] — 2F,+i [nl_Fjnl)^/^. As we show below, 

terms, particularly in 



{Fl 



-j+n"'i ->- -J L-J - '2Fj+i[n]Fj[n] 
this quantity is dominated by the F^ [ 
the higher-energy systems. 

A quick calculation shows that F^ [n] component is a sum 
of n variance terms and (n^ — n)/2 co-variance terms: 



Ff[n] 



J2 var(Af„) + 2 ^ 

i—j—n h^j—n 



E 

=/i+i 



cow{Mhj, Mi. 



(22) 
"We plot the square root of var(il/y) against j — i for 1/100 
filling in Figure [TO] Note once again the striking energy de- 
pendence. As a precursor to the Ff [n] data, the co-variance 
statistics are presented in Figure [^ as a function of i — h. 
Note the qualitative difference between the B = and B = 6 
cases in that the co-variance is negative without an applied 
magnetic field and positive with one. We have as yet no ex- 
planation for this difference. Another important feature is that 
the finite magnetic field co-variance is constant (though small) 
for nearly all i — ft,. Finally, the co-variance is noticeably larger 



for the lower-energy case. Aspects of the behavior of the co- 
variance of the Mij beyond RMT can, no doubt, be captured 
with the supersymmetric cr-model approach;^''*'* however, we 
are not aware of any results along these lines at this time. 

Figure E] depicts the resulting root-mean-square of Fj[n] 
as a function of n, averaged for all j in the energy band. All 
four systems considered in this paper are shown, and the quan- 
titative and qualitative differences between them are clear. 
Fj [n] becomes much larger for low filling cases, mainly due 
to the larger variance of Mij . At n = 400, Fj [n] in the B ^ 0, 
1/25 filling case is rising more rapidly with n than the zero 
field case because of the positive co-variance in the former. 

All four cases are noticeably different from the usual expec- 
tation for the behavior of Fj [n] . The expectation from RMT 
and random plane wave considerations is that Mij for differ- 
ent states nearby in energy are uncorrelated and have the same 
variance. Thus, var(Fj [n]) increases linearly at small n as un- 
correlated variables are added. However, the RMT modeling 
is only expected to apply up to SEij = Ei — Ej of order the 
Thouless energy i?Th — fiD/A. Beyond ETh, one can distin- 
guish in principle two energy ranges"^^ separated by the elastic 
scattering time r: (i) a first energy range i?Th < SEij < h/r 
corresponding to diffusive motion unaffected by the bound- 
aries, and (ii) h/r < SEij which is associated with the ballis- 
tic part of the dynamics. 

It is usually thought that the second of these energy ranges 
will be associated with the saturation of mesoscopic fluctu- 
ations. To understand the origin of this thinking, consider a 
quantity similar to Fj [n] but significantly simpler to analyze: 



i-i 



NAn]{r)^ J2 [\Mr)\'~l/A] 



(23) 



The magnitude of the fluctuations of Nj [n] (r) can be shown 
to be related to the probability of return of a trajectory to its 
original point r. However, the elastic time r sets a minimum 
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FIG. 12: (Color online) Standard deviation of Fj [n] for all four dif- 
fusive systems, averaged over all j in the energy window. The mag- 
nitude is surprisingly large. Predictions that Fj[n] would quickly 
saturate are not met; this is seen most clearly in the presence of a 
magnetic field where the increase is roughly linear even at n = 400. 



time before which no trajectory can return. As a consequence, 
no fluctuations are added by the energy range hJT < SEij. 

In the same way, it seems natural to expect that any meso- 
scopic fluctuation would saturate for energy larger than h/r. 
The systems we are investigating furthermore have been cho- 
sen in such a way that the elastic mean free path £ is not much 
smaller than the size L^ of the rectangle. As a consequence, 
the Thouless energy is not very different from the scale h/r. 

Thus, what we expect to see is a linear rise of the variance of 
Fj[n], followed by a saturation when SEij reaches an energy 
not much larger than the Thouless energy. This would also 
be the expected behavior for a ballistic system, provided one 
defines the Thouless energy as E^h = ^/tf where tf is the 
time of flight across system. For the parameters used here, the 
value of n at which we expect to see saturation, nxh, is 20-45 
for 1/25 fining and 10-25 for 1/100 filHng. 

What we observe, however, is a continued rise for all n, 
particularly with the inclusion of a magnetic field. In the case 
of 1/100 filling, we are near the bottom of the band, and the 
sum in Fj [n] provides a good estimate for summing over all 
the filled levels. Note that the continued linear increase of 
the standard deviation in the B ^0 cases requires correlation 
among the Mij . 

The only component not yet considered is the mean of 
Fj+i [n] ■ Fj [n], which one might expect to be about as large as 
F^ [n], thereby making var(i^j+i [n] — Fj [n] ) small. Although 
the statistic does similarly rise with n, it turns out to be consid- 
erably smaller than F^[n]: in both the 1/25 and 1/100 cases 
with magnetic field, this component reaches about one-fifth of 
the value of F^ [n] for the largest n. 

Finally, Figure^Jdepicts the quantity directly relevant for 
the Coulomb blockade peak spacing calculation, the root- 
mean-square of {Fj+i[n] — Fj[n]), averaged for all j in 
the spectrum. Showing remarkable similarity to the plot for 
a{Fj[n]) alone, including the order of magnitude. Figure [T3l 
contains odd features at the end of the 1/75 and 1/100 fill- 
ing plots that could either be noise caused by fewer eigen- 
functions being considered or real system effects involving the 
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FIG. 13: (Color online) The root-mean-square of (Fj+i [n] — Fj [n]) 
for all four diffusive systems, averaged over all j in the energy win- 
dow. This quantity is directly relevant to Coulomb blockade peak 
spacing. Although similar to those in Figure fT2l the curves are not 
simply \/2rms(Fj[7i]) because of correlations between Fj+i[n\ and 
Fj [n\ that are largest at low energy. As in Figure fT2l the magnitude 
is surprisingly large. 



very lowest eigenfunctions in the system. The most important 
feature of Figure^] however, is the sheer size of the Coulomb 
blockade-relevant statistic, even at high energies. When added 
to the peak spacing calculation, a statistic on the order of 1 
cannot help but cause major changes to system behavior. 

The end result of this section has been to demonstrate com- 
plex, unexpected behavior in the energy correlations of diffu- 
sive wavefunctions, behavior that we will demonstrate has a 
major effect on the statistics of the Coulomb blockade. 



V. COMPARISON TO ANALYTIC RESULTS 

We pause briefly from our main development to compare 
our results for the M^ to existing analytic results, particu- 
larly those from supersymmetric cr-model-'-^-^^ and random 
plane wave"*' calculations. The quantities we focus on are the 
mean and standard deviation of the IPR and off-diagonal Mij. 
These are given in the first six lines of Table I for our stan- 
dard four cases. Our interest, in particular, is in the deviation 
of these values from the universal values obtained in the sim- 
plest RMT - namely, that the mean should be integer (1, 2, or 
3 depending on the case) and the standard deviation should be 
0. We also give in the table the values for Fj [415] and the co- 
variance, the quantities showing the most unexpected results. 
No comparison of these will be made to analytic results, how- 
ever, because no such results exist. 

There are two caveats that one should bear in mind in mak- 
ing a comparison between the analytic results and our data. 
First, the analytic results are primarily for eigenfunctions that 
are close by in energy, within i^xii of each other There are 
no previous results, as far as we know, for our main finding 
that Fj {n) grows for large n because of coiTelations between 
eigenstates widely separated in energy. The second caveat is 
that the tr-model results are obtained in a somewhat different 
regime from our numerics: the cr-model assumes £ <C L while 
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Statistic 


1/25 Filling, B = 


1/25 Filling, B = 6 


1/75 Filling, B = 6 


1/100 Filling, B = 


= 6 


Mean IPR 


3.12 


2.14 


2.36 


2.29 


Mean Mi-i^i 


1.04 


1.05 


1.11 


1.09 


Mean Mij 


1.03 


1.03 


1.06 


1.05 


(7 (IPR) 


0.129 


0.082 


0.197 


0.203 


a{Mi-i,i) 


0.053 


0.039 


0.097 


0.092 


a{Mi,) 


0.044 


0.026 


0.056 


0.053 


rms(Fj[415]) 


0.67 


0.89 


N/A 


2.42 


rms(Fj+i [415] - Fj[415]) 


0.93 


1.10 


N/A 


2.07 


Mean cov{Mhj, Mij) 


-2.8 X 10"'' 


3.1 X 10"^ 


2.9 X 10"" 


2.8 X 10"^ 


Mean corr (Mhj , Mij ) 


-0.0016 


0.0055 


0.0097 


0.0106 



TABLE I: The energy correlation statistics most relevant to real mesoscopic systems are displayed for the four systems considered in this 
paper. In comparing system behavior, recall that the 1/100 filling system is at a lower W value (0.35 compared to 0.5) than the other systems, 
which was done to ensure it was diffusive. To see the effect of lowering energy while keeping the disorder value constant, one can compare 
1/25 and 1/75 filling. Note also the marked differences between the zero and non-zero magnetic field cases at constant energy. 



for our numerical results £ ^ L. 

The analytic approaches predict that three trends should be 
present in the data. First, the energy dependence: the devia- 
tions from simplest RMT should be proportional to the square 
root of the filling.^'"*' For our data, then, the deviations in 
the 1/100-filling column should be twice those at 1/25-filling 
(both at B — 6). We see that for the mean, this is certainly the 
case. For the standard deviation, the ratio ranges from 2.5 for 
the IPR to 2.0 for the far-off-diagonal case. Thus, this trend is 
quite reasonably obeyed by our data. 

Second, in the analytic approaches, the statistics of the IPR, 
Mii, is simply related to the statistics of the off-diagonal terms 
Mij . For the mean of the distribution, the cr-model approach 
yieldsi2i 



(IPR)b=o - 3 
{M,^j)b=o - 1 



3, 



(IPR) 



B^O 



-2 



(Mi^j)B^O - 1 



2 . 



(24) 



Our B~0 data is in good agreement with this result, for both 
the nearest neighbor values {i,i — 1) and the far-off-diagonal 
terms M^ . Note however that the ratio for our B ~ 6 data 
differs sharply from the above: for nearest neighbors the ratio 
is about 3 and increases to 5 or 6 for the far-off-diagonal terms. 
Turning to the standard deviation of the distribution, there are 
no (T-model results, but the RPW approach yieldsli 



rms(lPR)B=a 



s(M,^j) 



Ve , 



B=a 



rms(IPR)_B^o 

ims{M.,^j)B^o 



2 . 



(25) 



Good agreement is obtained for both of these ratios in the case 
of nearest neighbor terms, Afj i_i, but the far-off-diagonal 
terms show smaller fluctuations and so are not in agreement. 
Overall, then, the agreement between our data and this an- 
alytic trend is mixed: the results for the mean at zero field 
and the fluctuation of the nearest-neighbor terms is good, but 
those for the mean at non-zero field and the fluctuation of far- 
off-diagonal terms is poor. 

The third predicted trend is, of course, the relation between 
the results at zero magnetic field and those at non-zero field. 



For the mean, the d-model yields^ 



(IPR) 



B=0 



(IPR)s#o - 2 



= 3 



{Mi^j)B=o - 1 

{M,^j)B^a - 1 



= 2 



(26) 



In our data, however, the deviations in the mean are approx- 
imately independent of magnetic field. For the fluctuations, 
those of the IPR are expected to be larger by a factor of 3 in 
the presence of time-reversal symmetry.^' We see a factor of 
1.6 in our data. Thus there is a striking disagreement between 
the (T-model treatment and our numerics in terms of the effect 
of a weak magnetic field. 

To summarize our results with regard to the trends, some of 
the analytic predictions are seen in our data but others are not: 
the energy dependence checks, the behavior of the IPR com- 
pared to off-diagonal terms is mixed, and, most strikingly, the 
expected effects of breaking time-reversal symmetry are just 
not seen in the data. With regard to the latter, we emphasize 
that the simpler effects of breaking time-reversal symmetry, 
such as the change in level spacing distribution or distribution 
of |?/'(r)p, are certainly seen in our data, so the discrepancy 
here is not simply a matter of having applied too weak a field. 

To make a more exacting comparison of the data and ana- 
lytic results, we now compare the absolute magnitude of the 
deviation of the mean and the variance of the distribution. In 
order to do this, we need to first settle on a value for g, the 
dimensionless conductance, which is the main parameter con- 
trolling the expansion in the cr-model results. Standard ex- 
pressions exist for the conductance gdiff in the strongly dif- 
fusive Umit {£ ^ L) as well as for gbaih the ballistic con- 
ductance (i ^ L) assuming random scattering on the bound- 
aries. For our parameters we find at 1/25 filling gdiff == 36 
and 9baii = 16 while at 1/100 filling the values are g^is = 16 
and (^baii = 7.9. Our system is intermediate between these two 
limits. Surely the conductance cannot be larger than gbaii as 
this is the fundamental bound coming from the finiteness of 
the system. In fact, as a function of system size L while keep- 
ing other parameters constant, the conductance should cross 
over from gbaii oc L to the value gdiff which is independent of 
L. As our system is clearly in this cross-over region, use of 
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an interpolation formula appears necessary. We use the sim- 
plest such formula: g^^ — S^iff + ffbau- Thus the conductance 
values we use are gcS = H for 1/25 filling and ^cff — 5.3 for 
1/100 filling. 

With these values for gcS, the a-model expressions^ in the 
diffusive limit are easily evaluated. As the trends have been 
discussed above in detail, we give only a few representative 
values here. In the absence of time-reversal symmetry, we 
find 



(IPR) - 2 = 0.044 , rms(lPR) = 0.084 



at 1/25 filling and 



(IPR) - 2 = 0.096 , rms(IPR) = 0.18 



(27) 



(28) 



at 1/100 filling. Comparing with the Table, we see that for 
the fluctuations the numerical data are in good agreement with 
the diffusive cr-model results. In contrast, the deviation of the 
mean IPR from the RMT value is rather far off from the a- 
model predictions above. It is curious that once the scaling is 
done to change these values to those appropriate in the pres- 
ence of time-reversal symmetry [Eq. (I26lll. the fluctuations are 
not in agreement while the mean value of the IPR is. 

To summarize this section, we saw that there is already for 
the Mij statistics some substantial differences between our 
computational result and extensions of random matrix theory 
via the cr-model. Some of these differences are merely quan- 
titative, and might be explained by the fact that the regime we 
consider, where the mean free path is of order the system size, 
is not the one typically considered in cr-model calculations. 
Other differences, such as the existence of correlations among 
the Mij are qualitative, and thus less expected. 

The second kind of "integrated" statistics that we have con- 
sidered concerns the Fj's, which involve a further summation 
over states. For these quantities, we are not aware of any ana- 
lytic results. The fluctuations of the F/s show, however, (see 
Fig. ll3t an absence of saturation which is in total contradic- 
tion with intuition developed for simpler quantities. 

All these differences suggest that disordered quantum sys- 
tems may be a tougher nut to crack than previously thought. 



The peak spacing, with which we are most concerned here, 
is proportional to the second difference of the ground state 
energy with respect to electron number N: 

A'^En = E^,{N + 1) + Eg,{N - 1) - 2EgsiN), (29) 

which varies because of changing interaction effects as elec- 
trons are added and produces a peak spacing distribution. The 
simplest model of this distribution results from writing the 
ground state energy as the sum of the classical charging en- 
ergy and the energies of the occupied single particle states, 
known as the constant-interaction (CI) model. The ground 
state energy in this model can therefore be written in terms of 
the occupation numbers Uia- and one particle energies e,; as 






iNef/2C + J2^^'^^^ 



(30) 



where Uia- = 1 for the N lowest orbitals and zero otherwise. 
As a consequence, one gets the simple prediction 



eVCforoddiV, 



A^En 

A^En = eVC+(eAr/2+i-eAr/2)forevenAr. (31) 



The drastic odd/even difference is a quantum effect resulting 
from the spin of the electron. 

Although to our knowledge no quantum dot experiment has 
been conducted in the low-temperature, diffusive regime rele- 
vant to our study, all but one of the experimental results pro- 
duced so far for quantum dots show a marked disagreement 
with the prediction in Eq. ( 13 1> (for a recent discussion see Ref. 
47; the exception"*^ is the case of a quantum dot formed from 
a carbon nanotube). They show a wider peak spacing distri- 
bution and the lack of a pronounced odd/even effect. This has 
made it clear that the effect of the residual screened interaction 
between electrons is important in the description of Coulomb 
blockade experiments. A simplified but reasonably good ap- 
proximation for this residual interaction is a zero-range repul- 
sive potential 



l^s, 



-^^(r-r') 



(32) 



VI. APPLICATION TO THE COULOMB BLOCKADE 

A major way to probe the energy properties of electrons in 
a disordered media is by constraining groups of electrons in a 
quantum dot and studying the electrical transport of individ- 
ual electrons across the dot. The electrostatic charging energy 
being large, the dot is usually constrained to remain with a 
fixed number of electrons, which prevents current to flow at 
small bias voltage. This Coulomb blockade effect is essen- 
tially classical, and allows an applied gate voltage Vg to be 
adjusted so that the energy for N electrons is the same as that 
for A^ + 1, thereby inducing a finite conductance. The con- 
ductance through the dot as a function of Vg therefore forms a 
series of sharp peaks, the height and position of which encode 
information about the dot's ground statei^ 



(2) 

with Vq the total density of states (including spins) and Jg 
a parameter that can be taken as the Fermi liquid parameter 
/q ^ . We consider a case with moderate interactions: the value 
of Ts, the usual parameter to characterize the strength of inter- 
actions in an electron gas, is 1.5. For this value, Monte Carlo 



Ja) 



0.4, and we shall use this value of 



calculations'*' give /q" 
Js in what follows. 

It is known that in the absence of magnetic field, per- 
turbative calculations in this residual interaction should in- 
clude higher order terms (the so called Cooper channel) which 
makes the discussion significantly more involved.^'*'' We shall 
therefore restrict ourselves to the discussion of non-zero mag- 
netic field, for which a first order perturbation calculation is 
appropriate. In that case, the eigenstates are still Slater de- 
terminants characterized by occupation numbers riia = 0, 1. 
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FIG. 14: (Color online) Peak spacing distributions for the tliree sys- 
tems with 5 = 6 flux quanta and Js = 0.4. The integral under both 
the odd and even curve is normalized to 1, and both curves are shifted 
left a distance Js , as is customary. Note the wide spacing distribu- 
tion, the depression of the odd peak at Ja , and the lack of a strong 
odd/even effect, all increasingly apparent as energy is lowered. 



Noting that for a zero-range interaction like Eq. ( I32> the ex- 
change term exactly compensates the direct one for same spin 
electrons, their energy is given by Elni^a] = Eci[ni,a-] + 
Eri[ni,a] where i?ci["-j,o-] is the constant interaction contribu- 
tion Eq. ( I30t and 



-Eri [Ui 



J. 



2 S ^ 



^ijCr'^j^a' -^^^ij 



(33) 



is the residual interaction correction. For a given set of one 
particle energies e^ and wavefunctions t/i^ (thus fixing the 
Mij), the ground state is then obtained by minimizing E[ni,cr\ 
under the constraint that ^ n, ^r ~ N. Because of the resid- 
ual interaction term, this might not correspond to filling the N 
lowest one particle orbitals, and in particular may imply non- 
trivial (ie. not zero or one-half) spin S — ^ ■ ni^+ — J^i ^i-- 
(with + and — the minority and majority spins, respectively). 
If one assumes, however, the ground state occupation num- 
bers are the same as in the absence of interaction, Eq. J3U is 
just modified to 



A^E 



N — 



A'^En 



e^/C + Mn+i^n+i for odd iV, 

e /C + (eAr/2+l - eAr/2) 

+Fj+i [n] - Fj [n] for even N. 
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FIG. 15: (Color online) Spin distributions for the three systems with 
5 = 6 flux quanta and Js =0.4. Note that the average net spin 
increases as energy is lowered, and even at 1/75 filling, a spin of 1 is 
more likely than a spin of 0. The demonstrated spins are much larger 
than are predicted in most theories. 



Thus, the surprisingly large variance in both the Af^'s and 
Fj's may contribute greatly to the wider-than-expected peak 
spacing distribution observed^ experimentally. 

To confirm this, we applied the same eigenfunctions 
gleaned for the previous sections directly to a Coulomb block- 
ade calculation. All of the applied eigenfunctions have broken 
time-reversal symmetry and the interaction strength used is 
Js = 0.4. We show both the peak spacing distribution (Fig- 
ure [2} and net spin of the system (Figure [2} for all three of 
our relevant models. 

The CIh-RMT model, in comparison to Figure [2] shows 
a total concentration at zero for odd N and an asymmetric 
Wigner-Dyson-type distribution from about to 1 for even 
N. Other models, such as those relying on density functional 
theory, show departures from that basic structure, but typi- 
cally show a peak at Js for odd N and a comparatively thin 
distribution for even N. Our results predict that diffusive sys- 
tems should display a much wider peak spacing distribution 
and a disappearance of the odd/even effect. Both effects are 
increasingly marked at lower energies, which is also closer to 
the ideal inclusion of all energy states from the ground level 
up. 

The total spin is also much larger than most models predict, 
which is similarly likely due to interaction effects between 
distinct eigenfunctions. Thus, one prediction from this study 
is the presence of a large total spin in electrons constrained in 
a diffusive quantum dot at low temperature. 
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VII. CONCLUSION 

We have demonstrated that certain properties of disordered 
quantum dots are very different from expectations based on 
random matrix theory or random plane wave considerations. 
The key element in the properties which show these discrep- 
ancies is that they involve wavefunctions at different energies 
and an integral over space. They are properties which are both 
"ojf- diagonal" in energy and global in space. 

For these quantities, the mesoscopic fluctuations that we 
see are much larger than expected. Apparently the correlation 
among the wavefunctions amplifies the fluctuation effects. We 
emphasize that in making these statements we have been very 
careful to remain in the commonly defined diffusive regime. 
The mean free path (defined through the Born approximation) 
is less than the size of the system. And, with the exception 
of the width of the IPR distribution, all the properties which 
are either local in space or "diagonal" in energy are in good 
agreement with expectations. 

The explanation behind these unexpected results is largely 



open at this time. As briefly discussed in Section ITVAI one 
possible scenario would be that the observed statistics are the 
result of the presence of localized resonances. Further theoret- 
ical and numerical work would be needed to prove or disprove 
this suggestion, but if it held, it would give insight into the 
transition from the diffusive to localized regime in disordered 
quantum dots. 

In any event, the global off-diagonal quantities that we 
look at are exactly the quantities that are relevant to electron- 
electron interactions in finite systems. By looking at the addi- 
tion energy and ground state spin of our model quantum dots, 
quantities directly accessible to experiments in the Coulomb 
blockade regime, we showed that the unexpected statistics 
have a big effect on observable quantities. 
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